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Abstract. 

This  paper  presents  a general  overview  of  the  recent  works  dealing  with 
LES/DNS  carried  out  at  ONERA.  The  emphasis  is  put  on  developments 
aiming  at  extending  or  using  these  techniques  for  complex  flows  directly 
connected  with  practical  applications. 


1.  Introduction 

Unsteady  numerical  simulation  of  flows  of  theoretical  and  practical  interest 
is  more  and  more  popular.  It  allows  us  to  gain  a deeper  insight  in  involved 
physical  mechanisms,  and  practical  experience  shows  that  it  reduces  the 
constraint  on  the  turbulence  modeling  (when  compared  to  steady  RANS 
computations),  since  a significant  part  of  the  turbulent  motion  is  directly 
captured.  Main  drawbacks  are  that  the  solution  is  often  much  more  sen- 
sitive to  numerical  errors  than  steady  RANS  computations  and  that  huge 
amounts  of  grid  points  is  needed  to  obtain  a satisfactory  solution. 

The  rapid  increase  of  the  computing  power  makes  it  now  possible  to 
perform  DNS  and/or  LES  of  complex  flows,  which  cannot  be  considered  as 
” academic  test  cases”. 

In  this  paper,  we  will  discuss  some  applications  of  DNS  and/or  LES  to 
flows  of  practical  interest  performed  at  ONERA  : 

— Subsonic  flow  around  a two-dimensional  wing  near  stall. 

— Flow  around  a Delta  wing  with  vortex  breakdown. 

— Flow  around  a low-pressure  turbine  blade. 
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For  all  these  applications,  the  parameters  of  the  wind-tunnel  experi- 
ments (geometry,  Reynolds  number)  are  strictly  reproduced  in  the  simula- 
tions. All  these  flows  involve  a wide  range  of  physical  phenomena,  such  as 
transition,  separation,  boundary  layer,  wake,  etc.  As  a consequence,  they 
are  challenging  cases  for  LES  and  DNS,  since  they  cannot  be  treated  with 
one  of  these  approaches  only.  Because  of  the  high  value  of  the  Reynolds 
number  (up  to  several  millions  in  some  cases),  a full  DNS  cannot  be  em- 
ployed over  the  whole  computational  domain,  and  LES  must  be  used.  But 
transitional  region  must  be  treated  with  DNS  or  quasi-DNS  resolution. 
This  kind  of  complex  applications  thus  requires  the  development  of  hybrid 
approaches,  based  on  self-adaptive  subgrid  models  and  ad  hoc  numerical 
techniques,  which  must  be  able  to  provide  some  numerical  stabilization 
without  polluting  the  solution. 

New  developments  dealing  with  zonal  multiresolution  DNS/LES  tech- 
niques are  presented  in  Sec.  2.  Problems  dealing  with  the  control  of  numer- 
ical dissipation  are  discussed  in  Sec.  3.  Results  dealing  with  the  previously 
mentioned  applications  and  obtained  making  use  of  the  zonal  techniques 
are  presented  in  Sec.  4.  Then,  new  ways  to  reduce  the  computational  cost 
of  DNS/LES  are  discussed  in  Sec.  5. 

2.  Development  of  zonal  multiresolution  DNS/LES  techniques 

Most  of  the  known  works  dealing  with  LES  (see  [1]  for  a recent  review)  and 
DNS  on  structured  meshes  make  use  of  monodomain  numerical  algorithm, 
or  conformal  multidomain  approach.  This  results  in  a severe  limitation 
in  the  capability  to  concentrate  the  computational  nodes  in  areas  where 
important  physical  mechanisms  occur  while  keeping  the  total  number  of 
grid  points  as  low  as  possible.  This  limitation  renders  the  use  of  classical 
DNS/LES  approaches  very  expensive  for  complex  geometries,  where  these 
” important  areas”  are  very  small  when  compared  to  the  full  geometry.  This 
is  well  known  for  the  classical  example  of  the  turbulent  boundary  layer. 

In  order  to  alleviate  this  problem,  a zonal  multiresolution  approach 
has  been  developped  [18],  which  in  practice  appears  as  a non-conformal 
subdomain  technique.  From  the  theoretical  point  of  view,  this  approach 
leads  to  the  problem  of  the  coupling  of  two  solutions  with  different  filter 
lengths  at  the  subdomain  interface.  It  can  be  shown  that,  in  the  most 
general  case,  neither  the  instantaneous  field  nor  the  mean  or  rms  fields 
are  continuous  at  the  interface,  because  of  the  difference  in  the  filters.  As  a 
consequence,  classical  conservative  treatments  of  the  interface  are  no  longer 
valid. 

The  proposed  method  rely  on  the  definition  of  two  operators  at  the 
interface:  a restriction  operator  (from  fine  to  coarse  subdomains),  and  an 
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enrichment  operator  (from  coarse  to  fine  subdomains)  to  reconstruct  the 
high  frequencies  on  the  interface.  The  restriction  operator  is  defined  as  a 
filter  applied  to  the  high  frequency  solution.  The  enrichment  algorithm  is 
based  on  the  extrapolation  of  the  difference  between  the  solution  in  the 
two  subdomain  near  the  interface.  This  algorithm  was  demonstrated  to  be 
efficient  for  filter  length  ratio  as  high  as  4 for  wall  bounded  turbulent  flows 
(plane  channel,  2D  wings  or  delta  wing). 

The  use  of  this  technique  has  been  extended  to  the  definition  of  coupling 
between  2D  and  3D  subdomains  [15].  This  extension  makes  it  possible  to 
reduce  the  number  of  grid  points  for  flows  around  2D  profiles  by  using 
3D  subdomains  in  turbulent  zones  only  (boundary  layer,  wake),  while  2D 
domains  are  defined  everywhere  else. 


3.  About  the  control  of  numerical  dissipation 

Another  important  problem  which  occurs  for  high  Reynolds  number  com- 
plex flows  is  the  stabilization  of  the  simulation.  Daily  experience  shows  that 
subgrid  models  (even  the  most  dissipative  ones)  do  not  prevent  the  occu- 
rance  of  spurious  wiggles.  This  is  even  more  obvious  when  transonic  flows 
are  aimed  [7,  8,  9]. 

For  schock-free  flows,  a modified  version  of  the  AUSM+(P)  scheme  has 
been  developped  [15].  This  modification  is  based  on  the  introduction  of 
a wiggle  detector,  the  dissipative  part  of  the  AUSM  scheme  being  set  to 
zero  at  locations  where  no  wiggles  are  detected.  This  wiggle  detector  <f>  is 
computed  as  follows  in  cell  number  i : 

A!  _ f — 1 ^ (0i+ 2 — 4>i-\-l){(!)i+l  ~ 0i)  <0 

^ \ 1 otherwise  ' ' 


W*k  = 


l if  A«,t  + A»'<0  or  A^+A^co 
0 otherwise 


(2) 


The  sensor  is  then  defined  as  = Max(Wipk),k  = 1, .., 5,  where  ip  is 
the  primitive  variable  vector.  This  sensor  was  shown  to  yield  a very  large 
reduction  of  the  amount  of  numerical  dissipation,  preventing  a possible 
conflict  with  the  true  subgrid  dissipation  induced  by  the  models,  as  observed 
for  classical  upwind  schemes  [6,  17].  When  MILES  [2,  5]  computations  are 
wanted,  the  sensor  is  set  to  1 everywhere. 

When  shocks  are  present,  it  is  necessary  to  use  some  numerical  stabi- 
lization. High-order  ENO  filters  have  been  developped  [8],  which  lead  to  an 
important  cost  reduction  (up  to  a factor  15)  when  compared  to  classical 
ENO  schemes.  These  filters  are  localized  near  discontinuities  thanks  to  the 
use  of  a shock  detector  based  on  the  proposal  of  Ducros  et  al.  [4], 
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4.  Examples  of  application 

We  now  present  three  applications,  which  are  good  examples  of  recent 
achievements  of  LES/DNS  at  ONERA.  All  these  simulations  have  been 
carried  out  using  the  multidomain/multiresolution  technique  and  the  wig- 
gle detector  describe  above. 

4.1.  FLOW  AROUND  A TWO-DIMENSIONAL  WING  NEAR  STALL 

The  first  computational  case  is  related  to  the  flow  around  the  Aerospa- 
tiale A-airfoil  near  stall  (the  angle  of  attack  is  equal  to  13.3°).  The  chord 
Reynolds  number  is  equal  to  2.1  106.  This  case  was  selected  during  the 
LESFOIL  project,  and  details  of  ONERA’s  results  can  be  found  in  Ref. 
[15].  This  configuration  represents  a challenging  test  case  for  both  RANS 
and  LES  computations  because  of  the  complexity  of  the  suction  side  bound- 
ary layer:  the  adverse  pressure  gradient  induces  a laminar  separation  with 
turbulent  reattachment  near  the  leading  edge,  and  a massive  separation 
occurs  at  the  trailing  edge. 

The  best  results  have  been  obtained  when  the  transition  process  was 
directly  captured  by  the  simulation,  i.e.  when  the  grid  in  the  boundary 
layer  region  was  fine  enough  to  represent  the  small  scale  physics.  The  corre- 
sponding grid  contains  2048  points  on  the  suction  side,  yielding  A:r+  ss  100 
everywhere.  In  the  two  other  direction,  the  following  criteria  were  fulfilled: 
Ay+  < 20  and  A z+  < 2.  The  total  number  of  grid  points  is  7.2  millions. 
The  corresponding  number  of  grid  points  required  to  get  the  same  grid 
resolution  with  a monodomain  C mesh  is  21  millions. 

Good  results  have  been  obtained  using  the  Mixed  Scale  Model  [20,  11, 
12,  1].  Pressure  and  friction  coefficients  are  compared  with  experimental 
data  in  Fig.  1.  Computed  lift  and  drag  are  1.539  and  0.0250,  respectively. 
These  values  compare  favourably  with  experimental  results:  1.548-1.515 
and  0.0206-0.0308,  respectively. 

With  the  same  resolution,  MILES  results  are  less  accurate  but  still 
reliable,  mainly  because  the  transition  is  affected.  This  is  in  agreement 
with  the  results  of  Mossi  and  Sagaut  [17]  for  the  plane  channel  case. 

4.2.  FLOW  AROUND  A DELTA  WING  WITH  VORTEX  BREAKDOWN 

The  flow  around  a Delta  wing  at  high  angle  of  attack  with  vortex  breakdown 
is  an  example  of  3D  geometry.  The  Reynolds  number  based  on  the  chord 
and  the  angle  of  attack  are  equal  to  1.5  106  and  27°,  respectively.  We 
present  here  the  results  obtained  using  the  MILES  approach,  which  is  a 
good  candidate  for  this  kind  of  application  involving  massive  separation 
and  nearly-free  shear  layer  dynamics. 
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The  total  number  of  subdomains  is  equal  to  21,  corresponding  to  4 
millions  grid  points.  The  equivalent  number  of  grid  points  with  a conformal 
subdomain  technique  is  200  millions.  In  order  to  reduce  the  number  of  grid 
points,  a coarse  grid  was  used  in  the  near-wall  region  and  a wall-model  was 
employed  to  provide  boundary  conditions  on  solid  walls. 

A general  view  of  the  instantaneous  flow  is  shown  in  Fig.  2,  where  the 
vortex  breakdown  is  clearly  observed.  Mean  velocity  profile  is  compared 
with  experimental  data  in  Figs.  3 and  4,  showing  a satisfactory  agreement. 


4.3.  FLOW  AROUND  THE  T106  LOW-PRESSURE  TURBINE  BLADE 

The  last  example  deals  with  the  transitional  flow  around  the  T106  low- 
pressure  turbine  blade,  at  a chord  Reynolds  number  equal  to  1.6  105.  The 
inlet  flow  angle  and  the  exit  flow  angle  are  equal  to  37.7  and  -63.2  degrees, 
respectively.  The  pitch  to  chord  ratio  is  taken  equal  to  0.799.  The  compu- 
tational setup  corresponds  exactly  to  the  experimental  configuration.  This 
example  is  characterictic  of  applications  which  can  be  treated  using  both 
DNS  and  LES.  Extensive  DNS  simulations  have  been  performed  by  Durbin 
and  his  coworkers  [3]  on  the  same  geometry,  and  the  idea  was  here  to  reduce 
the  number  of  grid  points  while  preserving  the  quality  of  the  results. 

Once  again,  the  grid  was  fine  enough  to  capture  durectly  the  transition 
process:  Ax+  < 40,  A y+  < 20  and  A z+  < 1,  corresponding  to  the  use  of 
920  grid  points  on  the  profile.  The  total  number  of  grid  points  is  reduced  to 
1 million,  thanks  to  the  use  of  2D  subdomain  and  progressive  mesh  refine- 
ment in  the  3D  zone  along  the  profile  based  on  the  zonal  multiresolution 
technique.  The  total  number  of  subdomains  is  five  2D  domains  and  six  3D 
domains.  This  is  to  be  compared  to  the  25. 106  points  grid  used  for  DNS  at 
CITS. 

A general  view  of  the  instantaneous  flow  is  presented  in  Fig.  5.  The 
transition  process  and  the  topology  of  the  wake  are  clearly  identified.  Re- 
sults of  both  MILES  and  classical  LES  (using  the  Mixed  Scake  Model)  are 
compared  with  experimental  results  in  Fig.  6.  It  is  observed  that  a gen- 
eral very  good  agreement  is  recovered.  An  interesting  point  is  that  the  two 
approaches  yield  the  same  results.  This  can  be  explained  by  the  fact  that 
the  grid  resolution  is  very  good  (meshes  are  at  most  twice  larger  as  in  the 
DNS  grid),  leading  to  very  small  subgrid  effects.  The  discrepancy  with  the 
conclusion  drawn  from  the  A-profile  simulations  could  be  explained  by  the 
fact  that  the  size  of  the  mesh  in  streamwise  direction  is  twice  smaller  in 
the  blade  case  than  in  the  wing  computation,  and  the  Reynolds  number  in 
smaller. 
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5.  Further  reduction  of  the  computational  cost 

The  applications  presented  above  make  use  of  the  zonal  multiresolution 
technique,  but  they  all  rely  on  the  usual  LES/DNS  approaches.  Further 
reduction  of  the  CPU  cost  can  be  obtained  thanks  to  the  coupling  of 
DNS/LES  with  new  methods.  We  discuss  here  some  approaches  dealing 
with  the  coupling  of  LES  with  multilevel  algorithms  and  two  ways  of  per- 
forming RANS/LES  coupling,  which  have  not  yet  been  applied  to  complex 
configurations. 

5.1.  MULTILEVEL  ALGORITHMS 

A possible  way  to  reduce  CPU  cost  of  the  LES  approach  is  to  define  a multi- 
level algorithm  [13],  than  can  be  interpreted  as  a time-consistant  extension 
of  the  usual  multigrid  algorithm.  The  resulting  method  can  be  considered 
both  as  an  acceleration  technique  and  a new  way  to  close  the  filtered  Navier- 
Stokes  equations  without  making  use  of  the  subgrid  viscosity  approach.  In 
practice,  N levels  of  embedded  grids  are  defined,  each  one  corresponding  to 
a different  filtering  level.  The  key  points  of  the  method  are:  (i)  the  definition 
of  a cycling  strategy  and  (ii)  the  closure  at  each  level. 

Terracol  et  al.  have  proposed  a V-cycle  algorithm  [21,  22,  23],  based  on 
a fixed  or  self-adaptive  evaluation  of  the  integration  time  on  each  grid,  to- 
gether with  an  extension  of  dynamic  mixed  model  or  a hierarchical  closure. 
The  tests  show  a reduction  of  the  CPU  time  by  a factor  up  to  5 (compared 
to  monolevel  LES  on  the  same  grid)  without  noticeable  loss  on  the  accuracy 
of  the  results  for  the  plane  channel  flow  and  the  plane  mixing  layer. 

5.2.  ZONAL  RANS/LES  COUPLING 

A first  way  to  couple  RANS  and  LES  is  to  define  a zonal  approach,  as 
it  is  implicitely  done  in  the  Detached  Eddy  Simulation  approach.  From  a 
theoretical  point  of  view,  this  way  of  coupling  RANS  and  LES  appears  as 
an  extension  of  the  zonal  multiresolution  LES  algorithm  discussed  above 
in  the  paper.  This  analysis  is  conducted  by  Quemere  et  al.  [21,  19].  The 
main  new  points  are:  (i)  RANS  and  LES  solutions  do  not  correspond  to  the 
same  effective  filter  (that  can  be  analyzed  in  the  same  way  as  LES  effective 
filters  [14]);  (ii)  RANS  and  LES  solutions  do  not  necessarily  have  the  same 
dimension  and  (iii)  the  number  of  unknown  is  not  the  same  if  turbulence 
models  based  on  transport  equations  are  used. 

An  interface  treatment  for  zero-equation  SGS  model  and  two-equations 
turbulence  model  has  been  assessed  on  the  plane  channel  configuration, 
with  reliable  results  on  the  mean  velocity  and  the  resolved  Reynolds  stresses. 
Examples  of  computed  resolved  Reynolds  stresses  for  the  plane  channel  con- 
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figurations  are  shown  in  Fig.  7.  Two  configurations  have  been  investigated: 
near-wall  RANS  subdomain  with  LES  domain  at  the  center  of  the  channel, 
and  near-wall  LES  computation  with  RANS  domain  at  the  center  of  the 
channel. 


5.3.  ANOTHER  RANS/LES  COUPLING:  THE  NLDE  APPROACH 

The  last  point  discussed  in  this  paper  deals  with  the  reconstruction  of 
turbulent  fluctuations  around  a mean  profile  obtained  using  the  RANS 
approach.  An  approach  has  recently  been  developped  by  Labourasse  et  al. 
[21,  10],  which  belongs  to  the  Non-Linear  Disturbance  Equations  family,  as 
defined  by  Morris  and  his  co-workers  [16]. 

The  idea  is  here  to  reconstruct  locally  the  turbulent  fluctuations  for 
particular  purposes  (aeroacoustics,  aero-optics)  using  a LES-like  equations 
for  the  fluctuations,  with  source  terms  corresponding  to  the  coupling  with 
the  mean  flow.  Assessment  in  the  fully  general  compressible  case,  including 
unsteady  mean  flow,  on  the  channel  flow  configuration  is  given  in  Ref. 
[10].  These  results  also  show  that  the  NLDE  approach  is  more  robust  than 
classical  LES.  Examples  of  reconstructed  resolved  Reynolds  stresses  are 
compared  to  classical  LES  solutions  in  Fig.  8. 
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Figure  1.  Flow  around  a 2D  airfoil  near  stall  - Skin  friction  (Left)  and  Pressure  (Right) 
coefficients.  Line:  LES,  Symbols:  experiments 


Figure  2.  Flow  around  a Delta  wing  with  vortex  breakdown  (instantaneous  vorticity 
and  pressure  iso-surfaces) 
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Figure  S.  Flow  around  a Delta  wing  with  vortex  breakdown  - Mean  streamwise  velocity 
(Left)  and  mean  spanwise  velocity  (Right)  profiles  downstream  the  breakdown.  Line: 
LES,  Symbol:  experiments 


Figure  4 ■ Flow  around  a Delta  wing  with  vortex  breakdown  - Mean  wall-normal  velocity 
(Left)  and  mean  velocity  modulus  (Right)  profiles  downstream  the  breakdown.  Line:  LES, 
Symbol:  experiments 
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Figure  6.  Flow  around  a low-pressure  turbine  blade.  Mean  velocity  profile  (Left)  and 
RMS  velocity  profile  (Right).  Line:  LES,  Symbols:  experiments 
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Figure  7.  Zonal  RANS/LES  coupling  - Channel  flow,  ReT  — 590  - Streamwise,  span- 
wise  and  wall-normal  rms  velocity  fluctuations.  Near-wall  LES  treatment  with  RANS 
subdomain  at  the  center  of  the  channel;  line:  classical  LES,  symbols:  coupled  simulation 
(Left)  and  near-wall  RANS  simulation  with  LES  domain  at  the  center  of  the  channel; 
symbols:  classical  LES,  lines:  coupled  simulation  (Right).  Vertical  lines  note  the  position 
of  the  interface. 


Figure  8.  RANS/LES  coupling  - NLDE  Approach  - Channel  flow,  ReT  = 395  - Stream- 
wise,  spanwise  and  wall-normal  rms  velocity  fluctuations  (Left)  and  shear  stress  (Right). 
Line:  NLDE,  Symbol:  classical  LES 


